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Abstract. Tiie particle production process is reviewed, through which cosmic inflation can 
produce a scale invariant superhorizon spectrum of perturbations of suitable fields starting from 
their quantum fluctuations. Afterwards, in the context of the inflationary paradigm, a number 
of mechanisms (e.g. curvaton, inhomogeneous reheating etc.) through which such perturbations 
can source the curvature perturbation in the Universe and explain the formation of structures 
such as galaxies are briefly described. Finally, the possibility that cosmic vector flelds also 
contribute to the curvature perturbation (e.g. through the vector curvaton mechanism) is 
considered and its distinct observational signatures are discussed, such as correlated statistical 
anisotropy in the spectrum and bispectrum of the curvature perturbation. 



1. Introduction 

The standard model of cosmology at present is comprised by Hot Big Bang Cosmology and 
Cosmic Inflation. The cosmology of the Hot Big Bang accounts for the Hubble expansion of 
the Universe, the observed Cosmic Microwave Background (CMB) radiation, the primordial 
abundance of light elements (formed in the Early Universe through the process of Big Bang 
Nucleosynthesis) and the age of the Universe, which agrees well with the astrophysical estimates 
of the ages of the oldest globular clusters. In turn inflation overcomes or at least ameliorates 
certain fine-tuning problems regarding the initial conditions of the Hot Big Bang, namely the 
so-called horizon and flatness problems [1] [2]. 

What is cosmic inflation? In a nutshell, inflation is a brief period of superluminal expansion 
of space in the Early Universe. What inflation does is that it makes the observable Universe 
large and uniform. However, if the Universe were perfectly uniform there would be no structures 
like galaxies or galactic clusters, no stars with planets orbiting around them, no ... us. It is 
imperative, therefore, that there is a deviation from perfect uniformity, which can give rise 
to these structures. Indeed, we need a Primordial Density Perturbation (PDP) for structure 
formation to occur. It so happens that evidence for such a PDP exists as the latter reflects 
itself on the CMB through the so-called Sachs- Wolfe effect [3] which states that CMB light 
is redshifted when crossing growing overdensities. This effect directly connects the fractional 
amplitude of the PDP with the fractional perturbation of the temperature of the CMB: 
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Even though the PDP appears to be very small, numerical simulations of structure formation 



have shown that it is enough to account for the observed structure in the Universe. What is the 
origin of this PDP? Weh, it turns out that this too can be accounted for by Cosmic Inflation. 



2. Particle production during cosmic inflation 

To have an idea of how Cosmic Inflation is achieved consider the so-called Friedman equation, 
which is the temporal component of the Einstein equations for a homogeneous and isotropic 



spacetime. In flat space this equation reeds = Ip/m'p, where 
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Planck mas^ and H = a/a is the Hubble parameter corresponding to the rate of the expansion 
of the Universe, with a = a{t) being the scale factor of the Universe, parameterising the Universe 
expansion, and the dot denotes derivative with respect to the cosmic time t. Suppose now that, at 
some period in its early history, the Universe was dominated by an effective cosmological constant 
Agff. Then, /? ~ mpAes = constant, which means that the Hubble parameter H = Acs/3 is 
constant and, therefore, a oc e^*, i.e. space expands exponentially in time. Thus, inflation 
occurs when the Universe is dominated by an effective vacuum density. When inflation ends 
this density has to be transferred in the density of the thermal bath of the Hot Big Bang. This, 
in effect, amounts to a change of vacuurr0. Therefore, vacuum states during inflation are not 
necessarily vacuum states after inflation, but instead they can become populated. 

This process is called Particle Production and it is similar to the production of particles (in 
the form of Hawking radiation) on the event horizon of Black Holes [4J. Indeed, the cosmological 
horizon during inflation is an event horizon (of size ~ H^^) and can be viewed as an "inverted" 
(i.e. inside-out) black hole in the sense that nothing can escape being "sucked out" by the 
superluminal expansion. Virtual particle pairs, corresponding to quantum fluctuations, are 
broken up by the expansion and are pulled away to superhorizon distances, where they can 
no more find eachother and annihilate, becoming thereby real particles, giving rise to classical 
perturbations of the corresponding fields [5]. The amplitude of these perturbations is determined 
by the Hawking temperature 6(j) ~ Th, which for de Sitter space is Th = H/2tt [6J. Let us review 
now the particle production process in a more rigorous manner. 

The standard paradigm considers a real, minimally coupled, scalar field </> of mass m with 
Lagrangian density 

C = U,<t>d>^^-\m'<p'. (2) 

Using the above one obtains the equation of motion of the scalar field. The field is expected to 
become homogenised by the infiationary expansion as any inhomogeneities are inflated away. In 
this case the equation of motion becomes 
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where (j) = (j){t). Deviations from homogeneity are introduced originally only from its vacuum 
fluctuations. To follow their evolution we perturb the field from its homogeneous value as 
4> = (j){t) + 5(j){x,t). Then we obtain the equation of motion of the Fourier components of the 

field perturbations 5(j)k{k) = J 54>{x)e~^^'^dx. This equation reads 
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where k 



^ G is Newton's gravitational constant and we consider natural units where c — h — ks = 1- 

^ from the false vacuum of inflation, corresponding to Acff , to the true vacuum of the present, with A ~ 0. 



The next step is to promote the perturbations of the field to quantum operators defined as 



(2^) 



where a{k) and a}{k) are creation and annihilation operators respectively and we consider 
canonical quantisation with [a{k), a) {k' )] = {2tt)^ 5 {k — k'). The mode functions 6ipk{k,t) 
satisfy the same equation of motion as the Fourier components of the field perturbations 
54>k{k,t) because this equation is linear. We can solve this equation considering that in the 
subhorizon limit (k/aH — t- +oo) the solution matches the so-called Bunch-Davies vacuum 
which corresponds to the quantum fluctuations of a free scalar field in Minkowski spacetim 
Thus, the boundary condition reads 
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where r = —1/aH is the conformal time, which factors out the expansion of the Universe. Using 
the above, the solution of Eq. (jl]) for the mode functions is 



6^, = a~^/^/^ ^ ^^iJ Mk/aH) - e^^^ J.,{k/aH)], (7) 
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where Jy denotes Bessel functions of the first kind and u = — {^Y- To investigate particle 
production we evaluate the above solution in the superhorizon limit {k/aH — )• 0). We find 

s^, „-3/^,/Z£i^^ ('^\\ (8) 

^ V H sin(7ri/) T{l-v) \ k J ^ ' 

Using the above we can calculate the power spectrum Vs<i, = ^^\5(p\'^ in the superhorizon 
limit. We find 

87r|r(l -z/)|-2 /i7\2 / k ^ 3-21. 



l-cos(2^z.) \2tt) \2aH) ' 

Thus, we see that the scale dependence of the field perturbations is only included in the last term 
of the above. Considering a light field, i.e. m<Ci^=>i^— )-|, we see that the scale dependent 
term is eliminated and the spectrum becomes scale invariant. Indeed, for a light field the above 
spectrum can be written as 

where t/ = and we considered \'q\ <C 1. Thus, a light minimally coupled scalar field obtains 

a scale invariant spectrum of perturbations when they exit the horizon. The typical value of the 
field perturbations is 5(j) ~ ^Vs^ — /7/27r, i.e. it is determined by the Hawking temperature as 
expected. 

What is the fate of these perturbations after horizon exit? The perturbations can now be 
treated classically and their evolution is determined by the equivalent of Eq. ([3]) 

{5(^)+m{5cp) + m^{5(t))=Q. (11) 

Locally spacetime is flat. The subhorizon limit is well within the radius of curvature of spacetime during 
inflation so curvature can be ignored. Similarly, in this limit the momentum of the virtual particles is much larger 
than their mass fc/a 3> m so that the field can be considered effectively massless. 
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Figure 1. Schematic log-log plot of the evolution of the perturbations of a scalar field (enveloped 
by red solid lines) which follow the Universe expansion compared with the cosmological horizon 
(solid green lines). During inflation the cosmological horizon is an event horizon of constant 
physical radius ~ H~^. The perturbations of the scalar field start off as quantum fluctuations at 
subhorizon size, which are stretched to superhorizon scales by the (quasi) exponential expansion. 
After inflation, the size of the superhorizon perturbations continues to grow following the 
expansion of the Universe (i.e. proportional to a{t)) but the cosmological horizon (which is 
now a particle horizon) grows faster, with the speed of light. As a result, the perturbations 
reenter the horizon some time after the end of inflation. In the graph two different perturbations 
are depicted. They exit the horizon at different times, corresponding to the horizontal dashed 
lines. As a result they are inflated to different sizes so one is much larger that the other 
when they reenter the horizon. The size of the perturbations at horizon reentry is depicted 
by the vertical dashed lines. However, both the perturbations have approximately the same 
amplitude, determined by the Hawking temperature during inflation ^(t> = ^ = Th, which 
remains approximately constant. The fact that the scalar field perturbations retain the same 
amplitude even though they attain different sizes during their superhorizon evolution is the 
reason behind the scale-invariance of the perturbation spectrum. 

Considering that H ~ constant the solution of the above for a light field is 



Thus, we see that the perturbations of a light scalar field freeze-out after horizon exit. This 
guarantees that scale invariance is maintained since, regardless of the time of horizon exit (i.e. 
regardless of how large the size of the perturbations becomes by inflation), the perturbations 
have the same amplitude because H ~ constant until the end of inflation (see Fig. [T|). 

How are such field perturbations related with the PDF? The PDF arises because of the 
generation of a corresponding curvature perturbation (, as is discussed below. If the curvature 
perturbation is due to the perturbations of a light scalar field then their power spectra are 
proportional, i.e. oc Vs(j,- This means that they have the same scale dependence. The latter 





can be parametrised in the form of a power-law: Vc^ oc /c"' , so that 

Assuming that C, is due to a hght scalar field one obtains (cf. Eq. (|iup ) 

n, = l + 2r? + 0(e), (14) 

where e = —H/H^ <^ 1 quantifies the deviation from pure de Sitter expansion which we have 
ignored so far. Since, during inflation |7/| ^ 1 for a light scalar field, we find that, if the curvature 
perturbation is due to this held then Ug ~ 1 and the scale dependence of is eliminated. 
Thus, the PDP in this case would be approximately scale-invariant. Indeed, the latest WMAP 
observations suggest [8] 

= 0.963 ± 0.012 , (15) 

which agrees with the predictions of inflation with the quantum fluctuations of a light scalar 
field as the source of the PDP. Note that the observations deviate from exact scale invariance 
{hs = 1) at 1-a. This reveals some dynamics during inflation and agrees with the expectations 
of realistic inflation models. 

3. The curvature perturbation 

In general relativity the curvature of spacetime and the energy density of its content are 
interchangeable quantities, through the Einstein equations, depending on the choice of foliation 
of spacetime. Therefore, for the curvature perturbation it is useful to deflne a quantity C which 
is independent of such foliation (gauge invariant). This can be written as p] 

C = -V' - ^- , (16) 
P 

where the flrst term on the right-hand-side is the curvature perturbation in uniform density 
slices of spacetime and the second term is the density perturbation in flat slices of spacetime. 
We call C the (gauge invariant) curvature perturbation from now on. 

One can also define the power spectrum of the curvature perturbation (^^(x)) = d{ln k)V(^{k) 
which is given by the two-point correlator as 



(C(fe)C(fe')) = {27^f6^'\k + k')^^Vdk). 
The latest observations of WMAP give [8] 



(17) 



Vc_{ko) = (4.94 ± 0.09) x 10-^ (18) 

where k^ = 0.002 Mpc^^ is the pivot scale. The corresponding density perturbation (at horizon 
re-entry) is given by 

-) = ^ Cls = (1-98 ± 0.04) X 10-5, (19) 

P / H 

which, as mentioned, is the observed measurement of the PDPlfl 



* The subscript 'LS' refers to the last-scattering surface, where the CMB is emitted. 



Another useful quantity is the so-called bispectrum of the curvature perturbation, defined 
by the three-point correlator of C as follows 

iCmik'Kik")) = (27r)2<5(3)(^ + p + k")B^ik,k' ,k") . (20) 

The bispectrum can be related with the power spectrum is the following manner 

Bc(.ki,k2, h) = -^/nl[Pc(^i)^c(^2) + Pc(A:2)Pc(^3) + Pdk3)Pdh)], (21) 

where P(;{k) = ^V^k). 

The bispectrum is useful because it is exactly zero if the curvature perturbation is Gaussian, 
i.e. if it obeys Gaussian statistics. One expects that the perturbations of the fields which 
are generated during inflation from their quantum fluctuations are indeed Gaussian since the 
quantum fluctuations are random. Any non-Gaussian signal therefore would arise by the process 
which translates the perturbations of such fields, e.g. 5(j)^ into the curvature perturbation This 
process may be highly non-linear in which case Q would feature non-Gaussian statistics. The 
level of such non-Gaussianity is quantified in the bispectrum of Q by the so-called non-linearity 
parameter /nl- The latter also depends on the configuration of the three A;- vectors used to 
determine the bispectrum. 

The two most popular configurations used to determine /nl are 1) the equilateral 
configuration, where ki = k2 = k^ and 2) the squeezed configuration, where ki = k2 ^ k^. The 
WMAP findings for the values of /nl in these configurations are [S] 

/^l' = 26±140 and /^^^ = 32 ± 21 . (22) 

Note that there is a hint (at 1-a) of non-Gaussianity in the squeezed configuration. It is likely 
that non-Gaussianity may be detected in the near future by the observations of the Planck 
satellite, which are expected to improve the precision of the measurement of /nl by about an 
order of magnitude. It is important to stress here that the PDF is highly Gaussian. Indeed, 
remembering that V^{ko) ~ 10~^ and also i?^ cx ff^h'P'^ we see that, even if the observational 
upper bounds on /nl are saturated, non-Gaussianity in the PDP is tiny. This also agrees with 
the conjecture that the PDP is due to perturbations of suitable fields (e.g. light scalar fields) 
arising from their quantum fiuctuations during a period of infiation. 

4. The inflationary paradigm 

Before discussing the mechanisms through which the quantum fluctuations of suitable flelds 
can generate the curvature perturbation from inflation, it is necessary to briefly present the 
so-called inflationary paradigm. This is the typical manner in which inflation is modelled in 
particle cosmology. According to the inflationary paradigm the Universe undergoes inflation 
when dominated by the potential density of a scalar field, which is called the inflaton field. 

We return to the Lagrangian of Eq. ([2]) but this time instead of the mass term we consider 
a generic function V{4>), which corresponds to the potential density of the scalar field. In this 
case, Eq. 1^ has the form 

4> + 3H<j) + V'i(l)) = 0, (23) 

where the prime denotes derivative with respect to the scalar field (f). The above equation is 
similar to the equation of motion of a body sliding down the potential V and subject to a 
friction term determined by the rate of the Universe expansion H. Based on this analogy we 
will consider the field as rolling down the potential in field space. 
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Figure 2. Schematic representation of the inflaton potential V{(j)). During inflation the inflaton 
field 4> (represented by a circle) rolls down its flat potential in fleld space until it reaches a critical 
value (pend (corresponding to Vend = ^('?^cnd)) when slow-roll is terminated and inflation ends. 
Afterwards the fleld undergoes rapid oscillations around its VEV, corresponding to the minimum 
of the potential. The field decays while oscillating reheating thereby the Universe. 

Now, according to the inflationary paradigm, for inflation to take place the Universe must be 
dominated by the potential density of the inflation fleld, i.e. the kinetic density Pkinl*/*) of the 
fleld needs to be subdominant: 

VicI)) » Pkin(0) = ^0'. (24) 

This means that the field hardly moves and, therefore, its potential density remains roughly 
constant, providing thereby the effective cosmological constant Aeg = V/nip needed for 
infiationary expansion. When the condition p^in <C U is valid it is said that the field undergoes 
slow-roll and its variation (motion in field space) is overdamped by the friction term in Eq. ()23p . 
In this case one can ignore the acceleration term in Eq. (j23p which is then recast as 

3H<j) ~ -U'(0) . (25) 

Because the roll of the field is overdamped, the value of (p is very small, which means that the 
slope of the potential needs to be very small as well, according to the slow-roll equation above. 
Hence, the inflaton field corresponds to a flat direction in fleld space. 

Inflation ends when the slow-roll condition is violated, i.e. at some critical value (pcnd of 
the inflaton fleld, the potential becomes steep and curved such that pkin ~ V. Afterwards, 
the inflaton condensate oscillates around its vacuum expectation value (VEV). Such coherent 
oscillations correspond to massive particles (inflatonsjl, which eventually decay into the standard 
model particles that comprise the thermal bath of the Hot Big Bang (see Fig. [2]). This process 
is called Reheating and, if it occurs in a perturbative manner, it is usually modelled by adding 
a phenomenological decay term in the fleld equation 

4> + 3H<j) + r<j) + V'i(l))=0, (26) 



* with zero momentum since the field is homogeneous 



where T stands for the decay rate of the inflaton field. Reheating occurs when the density of 
inflation is transferred to the newly formed thermal bath, whose temperature is T^eh ~ V^^P^ 
called the reheating temperature. 



5. Mechanisms for the formation of the curvature perturbation 

We are ready now to discuss some ways that the perturbations of a light scalar field can give rise 
to the curvature perturbation in the Universe, which corresponds to the PDP that is responsible 
for structure formation. 



5.1. The inflaton hypothesis 

The traditional mechanism through which the curvature perturbation is generated employs the 
inflaton field itself for the job. This, so-called infiaton hypothesis, simply assumes that the field 
responsible for the formation of the curvature perturbation is the same field which determines 
the dynamics of infiation, i.e. the infiaton field. 

As mentioned in the previous section, to undergo slow-roll the inflaton must correspond to 
a flat direction in field space and therefore is characterised by V" <^ H^, i.e. it is a light scalar 
field. This means that, through the particle production process, it obtains a scale-invariant 
superhorizon spectrum of perturbations. As a result, the perturbed inflaton fleld will reach 
the critical value (j)end which terminates inflation at slightly different times at different points in 
space. Thus, inflation will continue a little bit more in some locations than in others (see Fig. [3]). 
This is the cause of the generation of the curvature perturbation d^. Indeed, it can be shown that 
the latter is determined by the difference in logarithmic expansion between the uniform density 
and the spatially flat slices, i.e. C = i5(lna) |10j . 

The PDP in this case is given by 
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(27) 



where the asterisk denotes the epoch of horizon exit of the inflaton perturbations. Since in 
the inflaton hypothesis a single degree of freedom (the inflaton) determines both the dynamics 
of inflation and the PDP, the latter can be written in terms only of the characteristics of the 
inflaton potential. Indeed, Eq. ()27p can be recast as [llj 
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which again is to be evaluated at horizon exit. Observational constraints on the amplitude of 
the PDP (cf. Eq. (jl9|) ) suggest that, if e is not extremely small, the energy scale of inflation is 
comparable to that of grand unification, i.e. V^^^ ~ 10^^~^^ GeV. This seems a natural scale 
to introduce new physics but it turns out that, for single field models, it requires fine-tuning of 
model parameters. 

The spectral index in the infiaton hypothesis is [llj 

n, = 1 + 2r/<^ - 6e , (29) 

where r/^ = m?pV" /V corresponds to the curvature of the potential along the direction of the 
infiaton field. If V{(t)) = \m^^(tP' as in Eq. ([2]), then t]^ = |(^)^ as discussed after Eq. (fTI#l . 

Under the infiaton hypothesis the generated non-Gaussianity is expected to be negligible 
with /nl ^ 1, i-e. below the limit of observability. This means that if non-Gaussianity is indeed 
observed in the PDP then all the single field infiation models are going to be falsified. 

^ we considered also the Friedman equation with p c^V since we have potential domination during inflation. 




Figure 3. Schematic representation of the way that the curvature perturbation is generated 
in the inflaton hypothesis. The inflaton field perturbations vary its expectation value typically 
by 6(p = H/2tt. This means that, while the inflaton field (represented by a circle) rolls down its 
potential in field space, its value is different throughout coordinate space, represented by the 
x-labelled axis. Thus, in some locations the rolling inflaton is ahead compared to others, which 
means that, in these locations, it will reach the critical value (pend that terminates inflation, 
somewhat earlier. At a given time the inflaton value is represented by the wiggly line which 
crosses the (p = (pend line in several places. Therefore, the Universe inflates more in some places 
than in others, because, at the given time that corresponds to the wiggly line shown, in some 
locations inflation has ended while in others it still continues. 

The beauty and the curse of the inflaton hypothesis is that it relies on a single degree 
of freedom, namely the inflaton fleld, to account for all the problems that inflation aims to 
address, i.e. provide the period of accelerated expansion that deals with the horizon and flatness 
problems and also generate the appropriate spectrum of curvature perturbations, which agrees 
with observational constrains on its amplitude and spectral index. As a result, models of single 
field infiation are typically overconstrained and suffer from substantial fine-tuning. This is why 
alternative hypotheses have been put forward for the generation of the PDP from inflation. 

5.2. The curvaton hypothesis 

This hypothesis assumes that the fleld responsible for the formation of the curvature perturbation 
has nothing to do with the dynamics of inflation, i.e. it is other than the inflaton field. This 
scalar field is called curvaton a [12] (see also Refs. [I3] [2] [15] [16] ) . In order for the curvaton 
to play this role it needs to obtain a superhorizon spectrum of perturbations during infiation. 
Thus, it needs to be a light field during inflation so that it can undergo particle production. 

It is evident that by introducing another degree of freedom, the fine-tuning problems of 
infiation model-building are substantially improved [T7][TB]. However, it must be stressed that 
the curvaton is not necessarily a new ad hoc addition to the theory, which is introduced by 
hand. Indeed, the fact that the curvaton is not connected to the infiaton sector allows it to 
correspond to physics at energy scales much smaller than that of inflation (e.g. the TeV scale 
which is accessible by LHC). Consequently, there exist a number of candidates for the curvaton 
which correspond to realistic flelds already present in simple extensions of the Standard Model. 



Prominent examples are a right-handed neutrino |19j . an MSSM flat direction |20) |21j or the 
so-called orthogonal axion in supersymmetric realisations of the Peccei-Quinn symmetry [22] [23] . 
Under the curvaton hypothesis the curvature perturbation is given by |24j 



(30) 



where f^o- 



rio-, where 0,^ = {pa/p)dcc is the density parameter of the curvaton field at 



the time of its decay after inflatiorlll. In the above Ca is the curvature perturbation attributed 
to the curvaton field, which is determined by the fractional perturbation of the field itselfH 
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where we used that Sa^, = i7*/27r and we assumed that near its decay the curvaton density is 
Po- ~ m'^a'^ , where mg- <^ is the mass of the curvaton field. The spectral index of the PDP 
in this case is [24] 

n, = 1 + 2r/, - 2e , (32) 



where rj^ corresponds to the curvature of the potential along the curvaton direction. 

Since by definition the curvaton should not affect the dynamics of infiation, during inflation 
we expect p^ <^ p- Thus the density parameter of the curvaton is extremely small and its 
contribution to the overall curvature perturbation (cf. Eq. ([30]) ) is also small. For the curvaton to 
significantly contribute to C we need to consider the period after infiation, when its contribution 
to the density can increase. 

In the simplest case when V{a) ~ ^m^a'^ the equation of motion for the curvaton is of 
the same form as Eq. ([3]). Since during infiation the field is light and rriu <^ H, the curvaton 
is overdamped and remains frozen in some value a^,. After infiation, however, the Hubble 
parameter decreases in time H(t) oc As a result, there will be a moment when rrio- ~ H{t), at 
which time the curvaton condensate will unfreeze and will begin coherent oscillations around its 
VEV. These oscillations correspond to massive particles (curvatons) whose density is diluted as 



Pa oc 



a ^, which is less drastic that the density of the radiation backgrounq£| p^ 



oc a 



Thus, the 



oscillating curvaton has a chance to dominate (or nearly dominate) the Universe before its decay 
(see Fig. H]). Thus, at curvaton decay Oq- can be as large as unity and the curvaton contribution 
to C can be substantial. Consequently, the curvaton imposes its curvature perturbation onto the 
Universe at (or near) its domination. 

Because of its spectrum of perturbations 6a the amplitude of the curvaton oscillations is also 
perturbed, i.e. it is larger in some places than in others. This means that the density of the 
oscillating field is perturbed too and the field will dominate the Universe at different times at 
different locations. This is what generates the curvature perturbation (cf. Eq. pip ). 

In contrast to the inflaton hypothesis, under the curvaton hypothesis non-Gaussianity can be 
substantial. Indeed, in this case [23] 

5 

•^NL = -7— , (33) 

which can be large if the decay of the curvaton happens before domination when <^ 1. In 
fact, the WMAP observations in Eq. (p2]) set the lower bound Oo- > 0.01. 



^ Clcr is also denoted as r in much of the hterature. 

* in spatially flat hypersurfaces 
9 



This is the radiation that was generated by the decay of the inflaton field after the end of inflation. 
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Figure 4. Log-log plot of the evolution of the inflaton energy density which decays into radiation 

(purple line) at the end of inflation (denoted by 'end') and the curvaton energy density 
(green line) (prompt reheating is assumed). During inflation, the curvaton density is negligible. 
After inflation p^ oc a~^. In contrast, p^^ remains constant (the curvaton is frozen at some 
value <T*) until m ~ H{t), when the curvaton unfreezes and begins oscillating, after which time 
(denoted 'osc') p^ tx a~^. At some moment (denoted 'dom') the curvaton density dominates the 
Universe {Q^ = ^^o- = 1) until, some time later (denoted 'dec') when it decays into the thermal 
bath of the Hot Big Bang. The dashed slanted line depicts the possibility that the curvaton 
decays before domination ((i^- = ^ 1), when substantial non-Gaussianity can be generated. 



5.3. Other mechanisms 

There are numerous other ways to generate the PDF from a superhorizon spectrum of scalar field 
perturbations. This section briefly reviews two of them, namely the inhomogeneous reheating 
and the end of inflation mechanisms. As with the curvaton mechanism, these mechanisms assume 
that the contribution to the curvature perturbation from the inflaton field itself is negligible. 



5.3.1. Inhomogeneous Reheating This mechanism generates the PDF by assuming that the 
infiaton decay rate T is modulated by another scalar field a [25] [26) . This scalar field is the one 
which undergoes particle production during infiation and which obtains thereby a superhorizon 
spectrum of perturbations. 

In this scenario we can ignore the perturbations of the inflaton field. Thus, we can consider 
that inflation is terminated at the same time throughout spacJ^. The inflaton then begins its 
coherent oscillations around its VEV until it decays when T ~ H{t). However, since in this case 
r = T{a) and a is perturbed in space, the decay rate is different at different locations so that 
the Hot Big Bang begins at different times. This is what generates the curvature perturbation 
which corresponds to the PDF. Roughly speaking we have 

5p 6T 6a 



at least in the observable Universe. 



In this scenario the spectral index is again given by Eq. ()32p . Non-Gaussianity in this model is 
quantified as [27] 

T'T 



/nl = 5 - 1 j , (35) 

where now the prime denotes derivative with respect to a. If T has a power-law dependence on 
a then /nl = 0(1)7 which is marginally observable. 



5.3.2. End of inflation This mechanism applies to a particular type of inflation model, the so- 
called hybrid inflation. Therefore, before discussing the mechanism a brief summary of hybrid 
inflation is necessary. 

Hybrid inflation couples the inflaton fleld cj) to another scalar fleld ip, called the waterfall 
fleld, in such a way that inflation is terminated by a phase transition which sends the waterfall 
to its VEV [25]. Since inflation is likely to be at the scale of grand uniflcation, in many cases 
the waterfall fleld is assumed to be the Higgs fleld of a Grand Unifled Theory (GUT). Thus, 
inflation is terminated by the breaking of grand uniflcation, i.e. the GUT phase transition. 

The scalar potential for hybrid inflation has the following form 

where M is the GUT energy scale, A is the self-coupling of the waterfall fleld and g is the 
interaction coupling between the waterfall and the inflaton. The potential V(},{4>) is responsible 
for the slow-roll of the inflaton and its precise form is not relevant to this discussioii^^l. From 
the above potential it is evident that the effective mass-squared of the waterfall field is 

m\ = gcp^ - XM'^ + \<i)'^ . (37) 

The above implies that there is a critical value of the inflaton 




such that if (/) » 0c then > and -0—^0. In this case, Eq. ([36]) becomes V = \XM^ + V^. 
The constant contribution to the scalar potential provides the effective cosmological constant 
for (quasi)de Sitter inflation: Agff ~ XM^/ntp. If however, cp < (pc then < and a phase 
transition occurs which results in ip ^ M, which gives the inflaton a large mass (through the 
interaction term) and so </> — t- 0. Assuming 1/^(0) = 0, after the phase transition y — t- and 
inflation ends. 

The end of inflation mechanism for the production of the PDF introduces an extra coupling 
between the waterfall fleld and another scalar field a [31j[32j. The scalar potential now becomes 

V{cP, V) = ^A(V' - My + \g<t?^l,^ + V^{ct>) + hia^i,^ (39) 

where h parametrises the strength of the interaction between a and the waterfall field ^ip. With 
this addition the effective mass-squared of the waterfall field becomes 

ml = g(p^ - XM^ + Xc/)"^ + ha^ (40) 



In supersymmetric versions of hybrid inflation the slow-roll potential is provided by radiative corrections and 
it is of the form V^oclnd) [^[50]. 



and the critical value which triggers the phase transition that ends inflation is now 

M<^) = i-M' - -aA , (41) 
\9 9 J 

i.e. it is modulated by the value of a. Therefore, if a is light during inflation, then it 
undergoes particle production and obtains a superhorizon spectrum of perturbations with typical 
magnitude 6a = H/2Tr. This means that the value of (j)c is also perturbed. Hence, the phase 
transition which terminates inflation occurs earlier in some parts of the Universe than in others 
depending on when the inflaton reaches the critical value 4>c- Consequently, the Universe inflates 
more in some locations than in others and this generates the curvature perturbation. As in the 
curvaton scenario, the spectral index in this case is given by Eq. ()32p |33] . 

6. Cosmic vector fields and the curvature perturbation 

Tantalising evidence exists of a preferred direction in the CMB temperature perturbations. In 
particular, the low multipoles in the CMB appear to be aligned |34j [35j at a level which is 
statistically rather improbable [3U|[37j. A preferred direction in the CMB cannot be accounted 
for if inflation and the generation of the PDF is due to scalar fields only. Moreover, despite their 
abundance in theories beyond the Standard Model, scalar fields have not been observed as yet. 
If the Higgs field is not found in the LHC, the credibility of the ubiquitous usage of scalar fields 
in cosmology will be shaken. 

Until recently only scalar fields have been considered both as responsible for the dynamics 
of inflation and also for the generation of the observed PDF. However, in the pioneering work 
of Ref. [38j the possibility that a vector field contributes in the generation of the PDF was first 
considered. Since then, a number of attempts have been made to investigate the role and the 
implications of cosmic vector fields in the generation of the PDF (see Ref. |39j and references 
therein) . 

Why vector fields were not considered originally for the generation of the PDF? After all, 
these are fields which are similar to the massive gauge bosons, which have indeed been observed 
in LEP. The main difficulty had to do with the inherent anisotropic nature of vector fields. 
Inflation would homogenise a vector field condensate, and a homogeneous vector field picks 
up a preferred direction in space. Thus, if this vector field condensate were to dominate the 
density of the Universe (so that it can affect the expansion and generate the PDF) it was 
felt that the anisotropic stress generated would lead to strongly anisotropic expansion which 
would be impossible to reconcile with the predominant isotropy of the CMB. However, as is 
discussed below, this difficulty can be circumvented in the vector field remains subdominant 
during inflation. Still, if (in analogy to scalar fields) light vector fields were needed for a 
scale invariant spectrum then a more subtle problem arises. Massless Abelian vector fields 
are conformally invariant and they cannot undergo particle production during inflatiorF^. Thus, 
for light vector fields particle production is expected to be suppressed. Nevertheless, there are 
numerous mechanisms which break the conformality of vector fields so that particle production 
can be efficient. These mechanisms are model dependent which means that they can have 
distinct observational signatures as is discussed below. 

7. Particle production of vector fields during infiation 

Suppose that there is a suitable theory which breaks the vector field conformality during inflation. 
How are we to study particle production of the vector field? We follow the same recipe as in 
Sec. [2] assuming that the inflationary expansion remains isotropic. 



They view the Universe expansion as a conformal transformation to which they are insensitive. 



Firstly we perturb the vector field around its homogenised value Afj_ = ^^(t) + 6Afj_{x,t). 

Then we Fourier transform the perturbations 6Ai{k) = J 6Ai{x)e~^^'^dx and we obtain the 
equations of motion for the Fourier components in the theory which we are considering^. Then 
we promote the perturbations of the vector field into quantum operators 

6Mx,t) = J ^Y.[^f^>^(k)6Ak{k,t)e''^-^ + e^*a^^^^^ , (42) 

where 6Ak are the mode functions and we consider again canonical quantisation 
[ax{k),a\,{k')] = (27r)^5(A: — k')6\x'. In the above we have introduced the polarisation vectors 

e^ = -^(l,i,0) e^ = -^(l,-i,0) ell = (0,0,1), (43) 

where L,R denote the left and right transverse polarisations and || denotes the longitudinal 
polarisatiorF^. 

As with the scalar field case, the mode functions of the above expansion are expected to 
satisfy the equations of motion for the Fourier components of the vector field perturbations. 
The form of these equations depends on the theory which breaks the conformality of the vector 
field. To solve them we again consider vacuum boundary conditions in the subhorizon limit 
(k/aH — )• oo). They are (41j : 

ik/aH ik/aH 

5Ab^ > and bA > 7-^. (44) 

We see that the boundary condition is again the Bunch-Davies vacuum but the longitudinal 
component it is multiplied by the Lorentz boost factor 7 = ^ = ^ (^)^ + which takes us 

from the frame with k = where all components are equivalent, to the frame of momentum fcP^ 
After solving the equations of motion we evaluate the solutions in the superhorizon 
limit {k/aH 0). Then we can obtain the power spectrum for each polarisation as 
V\{k) = ^{SA^l"^, where 6A^ is evaluated in the superhorizon limit. Because we have three 
degrees of freedom there are three possibilities: 

Case A: V\\ ^Vl^Vr^ V\\ 
Case B:VL = VRi- V\\ 
Case C: Vl=Vr = V\\ 

Even though Case A appears to be the most generic, in practice it is Case B that is the 
most common because usually the theories which break the vector field conformality are parity 
conserving. We have isotropic particle production only in the special Case C. 

What would happen if the spectrum of the produced perturbations of the vector field affected 
the curvature perturbation? Since particle production is in general anisotropic we expect an 
anisotropic contribution to C,. This contribution can be parametrised as follows. For the power 
spectrum we can write 



V^{k)=V'^^{k)[l + g{A-kf...\, (45) 

We focus on the spatial components of the cosmic vector field, as they are the ones which would generate 
anisotropy. Also, the temporal component of a homogeneous massive Abelian vector field is zero. 
The index A in Eq, (|42|l runs over these values. 

Note that for a massless vector field the longitudinal component is unphysical. 




Figure 5. Patterns in the CMB temperature perturbations 
which may arise from statistical anisotropy in the CMB 
spectrum. The left panel shows an isotropic signal, 
while the middle and right panels show patterns due 
to statistical anisotropy along the vertical and horizontal 
direction respectively. The figure is taken from Ref. [44j . 



where A = A/A and k = k/k with A = \A\ and the ellipsis denotes higher order contributions. 
The anisotropic part of the spectrum corresponds to a new observable, namely Statistical 
Anisotropy, which amounts to direction dependent patterns in the CMB [Jl] (see Fig. [5]). 
Statistical anisotropy is quantified by the so-called anisotropy parameter g in Eq. ()15|) . The 
WMAP observations set a surprisingly weak upper bound on the anisotropy parameter: g < 0.3 
. Thus we see that, at present, observations allow up to 30% statistical anisotropy in the 
spectrum of the PDF. The forthcoming data of the Planck satellite will improve precision by 
an order of magnitude and are likely to detect statistical anisotropy. If this will be so then a 
cosmic vector field will have to be involved in the generation of the PDP during infiation. 

In a similar manner one can parametrise the statistically anisotropic contribution from a 
vector field in the bispectrum. Indeed, we can write [H] 06] 

fNL = f^l{l + GAl---), (46) 

where /j|£ denotes the isotropic part and A_i_ is the projection of A onto the plane of the 
three /c-vectors which are used to define the bispectrum. Note that, as non-Gaussianity has 
not been observed yet, there are no bounds on ^. If ^ > 1 this means that non-Gaussianity is 
predominantly anisotropic (even though is not). If non-Gaussianity is indeed observed (e.g. 
by the Planck satellite) and no angular modulation of /nl is found on the microwave sky then 
all models which predict ^ > 1 will be falsified. It is important also to stress that the directions 
of statistical anisotropy in the spectrum and the bispectrum are correlated (determined by the 
direction of A), which is a smoking gun for the contribution of a vector field to the PDP [35]. 

In the Cases A and B the spectra of the superhorizon perturbations of the vector field are 
predominantly anisotropic. Since is isotropic at least at the level of 70% the contribution 
of the vector field to the PDP in these two cases has to be subdominant and its significance is 
only the possibly observable statistical anisotropy in the spectrum and bispectrum. Thus, in 

In Ref. 1^ statistical anisotropy with g = 0.29 ± 0.03 is claimed to be detected at a level of 9-a but the authors 
acknowledge that, because the direction is suspiciously near the ecliptic plane, this is probably due to an unknown 
systematic error. Hence we treat this number as an upper bound. 




the Cases A and B we still require some other, isotropic source (such as a scalar field) to provide 
the dominant contribution to the PDP. On the other hand, in Case C the perturbation spectra 
of the vector field are isotropic. In this case, the vector field alone can generate C and no input 
from any scalar field is necessary. Note that, in this case, the vector field may also produce 
statistical anisotropy if the Vx are not exactly the same but differ slightly, albeit by no more 
that 30%. 



8. Models for particle production of cosmic vector fields 

In this section we discuss a couple of proposals for the generation of a flat superhorizon spectrum 
of vector field perturbations. 



8.1. Non-minimal coupling to gravity 
Consider the following theory 

C = --^F^uF^^" + \{m^ + aR)A^A^', (47) 

where -F^^ = dfj_Ay — dyA^ is the field strength of the Abelian vector field A^, m is its bare mass, 
a is constant and R is the Ricci scalar. From the above we see that the vector field has effective 
mass-squared m\ = m? + aR. 

Following the procedure outlined in the previous section, we obtain the following_solution for 
the mode functions of the transverse perturbations of the physical vector field [47 



[gS( 

B 



5^.^ = a-'I^J^^^^^[Uk/aH) - e^-'^J.^k/aH)]. (48) 



The above is identical to Eq. ^ with the crucial difference that i'=\J\ — (^7^)^- Using 
this solution we obtain the following expression for the power spectrum of the transverse 
perturbations in the superhorizon limit 

l'-Ls(2Ji) (2^) (2^) ' ^^^^ 



which again is identical to Eq. ([9]) but for the different value of v. From the above it is evident 

3 
2- 



that a scale invariant spectrum is attained for |. This translates into a requirement for 



both a and m as follows. 

In a spatially flat, homogeneous and isotropic spacetime, the scalar curvature is 
R = 3(3t« - l)H'^, where w is the barotropic parameter of the content of spacetime. During 
inflation, ~ —1 so that ~ — 12 . This means that we can attain a scale invariant transverse 
spectrum if 1) a = h and 2) m H, i.e. the vector fleld is a light field with a negative effective 
mass-squared ^Zlij! Under these conditions the transverse spectra are indeed scale-invariant 
and are given by Vr^l = (^)^> exactly as a massless scalar field. 

Employing the same conditions for a and m one obtains the following solution for the 
longitudinal component of the physical vector field 



5A 



1 

V2 



ik/aH 

(50) 



^"^ In a flat homogeneous and isotropic Universe the spatial components of the physical (in contrast to comoving) 
vector field are Ai/a. 

The non-minimal coupling to gravity a = ^ has a special property. It turns the conformally invariant massless 
Abelian vector field equivalent (in the sense that its field equations are the same) to a (set of two) minimally 
coupled massless scalar field(s). If applied to a massless scalar field it renders it conformally invariant. In this 
respect it appears to refiect a deeper symmetry and therefore it may be a natural value for a. 



Clearly the above is totally different from Eq. (|48p . In the superhorizon limit the power spectrum 
of the longitudinal component is Vw = 2(^)^ = This means that we are in Case B since 

^11 > 'Pl,r and Vl = Vr because the theory is parity invariant. Thus, particle production in this 
theory is anisotropic at a level of 100%. Therefore, the vector field cannot alone be responsible 
for the PDF. Its contribution to C, has to be subdominant and it can only be the source of 
statistical anisotropy. Before concluding, we should point out that this model has been criticised 
for suffering from instabilities (ghosts) [lS][12j (see however Ref. [511]). 

8.2. Varying kinetic function and mass 
Now, consider the theory 

>C = -J/F^.F'^- + ^m%A^ (51) 

where F^i, = d^Ay — d^A^, f = f{t) is the kinetic function and we also consider m = m{t) > 
during inflation. The Maxwell-type kinetic term in combination with the positive mass-squared 
guarantees the stability of the model [51j and therefore is motivated even if the vector field is 
not a gauge boson. Note also that a massive Abelian vector field is renormalisable even is it is 
not a gauge field [52] . 

The solutions for the mode functions of the field perturbations are too complicated to 
reproduce here. It suffices to say that scale invariance in the transverse components requires 
the kinetic function to scale with the expansion as / oc a^^^^ and also the physical vector field 
to be light at horizon exit M <^ H, where M = . Scale invariance for the longitudinal 

component additionally requires m (x a [39] [53] . 

If the vector field is a gauge boson then / is the gauge kinetic function, which is related 
to the gauge coupling as / ~ 1 /e^ . This means that only the case when / oc a~'^ is acceptable 
because only then the gauge field remains weakly coupled during inflation. Note that the model 
is naturally realisable in supergravity theories where / is a holomorphic function of the scalar 
fields of the theory f{4>i) [53]. In general, Kahler corrections to the scalar potential result in 
masses of order H for the scalar fields [55jij, which are therefore expected to be fast-rolling down 
the potential slopes during inflation, causing significant variation to /. Indeed, for a power-law 
dependence of the gauge kinetic function to the scalar fields, it is easy to show that f / f ^ H [39j, 
i.e. / has a power-law dependence on a as assumed in this model. Now, if / is modulated by the 
infiaton field itself, then it can be shown that, under fairly general conditions, the backreaction 
to the inflaton's roll renders the scaling / oc an attractor solution j56j. Similarly, for a gauge 
boson, m can be modulated by a fast-rolling Higgs field with tachyonic mass mn = 2i7 |39| . 

Given the conditions / oc a~^^^ and m cc a the power spectra for the transverse and 
longitudinal components depend on whether the vector field remains light until the end of 
infiation or not (note that M = M{t)). In particular we find [39j [53j 

M««: f. = fK=(^)' and f„ = (|)^ (f )^ (52) 
M»«: p, = p„ = ^„ = ^(ii)^(^)\ (53) 

From the above we see that, if the vector field remains light until the end of infiation, 
V\\ ^ Vl = Vr, i.e. we are in Case B. Therefore, particle production is strongly anisotropic 
and the vector field contribution to the FDF has to be subdominant but it can still give rise 
to substantial statistical anisotropy. In contrast, if the field becomes heavy by the end of 
infiation then particle production is isotropic and the vector field can be solely responsible for 

This is the source of the famous 77-problem of inflation, since slow-roll requires \ri\ — |(f )^ < 1 for the infiaton. 



the generation of the PDP. Note that, because we need the field to be Hght when the cosmological 
scales exit the horizon, for it to become heavy we need M > during inflation, which is possible 
only in the case when / cx a~^. If / oc then M = constant and we have M <^ H throughout 
inflation. But if / oc a~'^ then we have M tx so that we may end up having M ^ H at the 
end of inflation even though we started of with a light field at horizon exit. In the latter case 
the vector field begins coherent oscillations a few exponential expansions (e-folds) before the 
end of inflation. Assuming that at the end of inflation / — )• 1 and the field becomes canonically 
normalised, we find that there is ample parameter space for Case C to be realised, namely: 
1 < m/H < 10^ I32][S3]. 

9. The vector curvaton paradigm 

Through the examples in the previous section it is evident that a scale invariant spectrum 
of perturbations (isotropic or not) of a vector field can indeed be generated if one assumes 
some theory which appropriately breaks the conformality of the vector field. However, in order 
for these perturbations to affect or even generate the PDP the vector field needs to affect the 
Universe expansion, i.e. its density should become dominant (or nearly dominant) at some point. 
But, even if particle production is isotropic the homogeneous zero- mode condensate is not. How 
can we avoid excessive anisotropic stress when the vector field dominates the expansion? It 
turns out that this is possible if the vector field plays the role of the curvaton. 
Consider a minimally coupled massive Abelian vector field for which 

C = -\f^uF^'' + \m^A^A^', (54) 

where F^j^ = d^A^, — dyA^j^ and m = constant > 0. It is clear that both models in the previous 
section eventually approach the above theory, when m? ^ ii ~ and also after the end of 
inflation when / — >• 1 and m = constant. The vector field condensate is homogenised by inflation 
so that A^ = A^{t). In this case it can be shown that the temporal component is zero and the 
spatial components satisfy the following equation [38] 

+ HAi + m^Ai = , (55) 

which is similar to Eq. ^ for the scalar field. The energy- momentum tensor of this theory can 
be written in the form [38j 

= diag(/3A, -P±, -P±,+Pi.) , (56) 

where 

PA = Pkm + V Pkhi = -^F^uF'"' 

and (57) 
P± = Pkin -V V = -^m^Af.A''. 

Eq. (|56p is reminiscent of a perfect fluid with the crucial difference that the pressure in the 
longitudinal direction is of opposite sign compared to the transverse pressure. Thus, it seems 
that if the homogeneous vector field were to dominate the Universe it would indeed generate 
excessive anisotropic stress. Therefore, we assume that the vector field remains subdominant 
during infiation so that (quasi)de Sitter expansion is not spoilt. 

After inflation the Hubble parameter decreases until m > H{t). When this happens one can 
ignore the friction term in Eq. (j55|) . which therefore suggests that the vector field condensate 
starts rapid (quasi) harmonic oscillations. It is easy to show that, during a Hubble time, on 
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Figure 6. Log- log plot of the evolution of the inflaton energy density which decays into radiation 
(purple line) at the end of inflation (denoted by 'end') and the vector curvaton energy density 
PA (green line) (prompt reheating is assumed). During inflation, the vector curvaton density is 
negligible (its evolution depends on the model which breaks its conformality) . After inflation 
p^ oc a~^. Similarly, pA decreases also as radiation after inflation when the vector field is light 
(in contrast to the scalar curvaton case where p^ remains constant when the curvaton is light, cf. 
Fig. However, when m ~ H(t), the vector curvaton becomes heavy and begins oscillating, 
after which time (denoted 'osc') pA oc a^^- At some moment (denoted 'dom') the vector curvaton 
density dominates the Universe until, some time later (denoted 'dec') when it decays into the 
thermal bath of the Hot Big Bang. As the vector curvaton is rapidly oscillating it does not give 
rise to anisotropic stress at domination. The dashed slanted line depicts the possibility that the 
vector curvaton decays before domination (Qa ^ 1)) when substantial non-Gaussianity can be 
generated. 

average pkki ~ which means that, over many oscillations, the average pressure is zero: p± = 
[38]. Thus the oscillating vector field condensate behaves as pressureless, isotropic matter. Its 
density, therefore, scales as pA oc which is not as drastic as the radiation background, which 
scales as p^ oc a""^. Therefore, the oscillating vector field can gradually increase its density 
parameter and come to dominate (or nearly dominate) before its decay (see Fig. [6]). Because it 
is isotropic, it dominates without causing any excessive anisotropic stress, so that the Universe 
expansion remains isotropic |38j . 

From the above we see that a massive Abelian vector field can follow the curvaton scenario 
and play the role of the curvaton without problem. Indeed, the perturbations of the vector field 
imply a perturbation in its local density pA {x) , which means that in some locations it dominates 
the Universe earlier than in others. This is how it can generate the curvature perturbation 
according to the curvaton mechanism. It is important to note that, since pA is a scalar quantity, 
the curvature perturbation generated is scalar and not vector in nature. Also, note that the 
perturbations of the vector field 6Ai follow a similar equation to Eq. (j55|) . which means that 
they too are rapidly oscillating and do not introduce anisotropic stress. Thus, at domination the 
perturbations do not cause anisotropic expansion either, not even in a small scale. If particle 



production of the vector field is anisotropic then there are direction dependent patterns in the 
amplitude of the oscillating zero mode that lead to statistical anisotropy in the PDP [41] . 

In Ref. statistical anisotropy in the bispectrum in the vector curvaton model was 
investigated. It was found that it manifests itself only quadratically, i.e. the expansion in 
Eq. (jlSl) is truncated to /nl = Ml (l + ^^l) ■ ^ 

was evaluated in the two models discussed in 

the previous section. In the non-minimally coupled to gravity model Q'^"'^ = | in the equilateral 
configuration, whereas Q^"^^ = 1 in the squeezed configuration. This means that the angular 
modulation of /nl is prominent and should be detected if non-Gaussianity is found. In the 
varying kinetic function and mass model Q""^^ = > ^""^^ = (ig-)^ > 1 ^^en M < H, 

i.e. if the field remains light until the end of inflation. Here we see that non-Gaussianity is 
predominantly anisotropic. Should no angular modulation of /nl be observed this possibility 
will be excluded. However, if M ^ i.e. if the field becomes massive by the end of inflation, 
particle production is isotropic (Case C) and /nl = /nl ~ 5/40^, which is identical to the scalar 
curvaton case (cf. Eq. (f33j) )F^ 

The vector curvaton is an elegant mechanism for vector fields to contribute of even generate 
the PDP since it does not need to couple the fields to the inflaton sector. However, it is by no 
means the only way that a vector field can affect the curvature perturbation. For example, in 
Ref. [57], the end of inflation mechanism is employed (see Sec. 15.3.2]) . where, instead of coupling 
the waterfall field tp of hybrid infiation to a scalar field cr, the authors considered introducing 
a coupling of the form /SV = ^hA^A^i/^"^ . This model too produces distinct observational 
signatures. For example, in this model /nl = /nl(-'- + + G'A'^), where the maximum value 
for Q' in the squeezed configuration is independent from the model parameters: G'y^^^ ~ j- 

10. Conclusions 

Cosmic structure originates from the growth of quantum fiuctuations during a period of cosmic 
inflation in the Early Universe. The particle production process generates an almost scale 
invariant spectrum of superhorizon perturbations of suitable fields, e.g. light scalar fields. 
These perturbations give rise to the primordial density /curvature perturbation via a multitude 
of mechanisms (e.g. inflaton, curvaton, inhomogeneous reheating etc.). Observables such as the 
spectral index Ug or the non-linearity parameter /nl will soon exclude whole classes of models. 
Indeed, the Planck satellite observations are expected to increase precision up to /nl = 0(1)- 

Recently, the possibility that cosmic vector fields contribute or even generate the curvature 
perturbation (e.g. through the vector curvaton mechanism) is being explored. If it is so 
then vector fields can produce distinct signatures such as correlated statistical anisotropy in the 
spectrum and bispectrum of WMAP observations allow up to 30% statistical anisotropy in 
the spectrum but the Planck satellite mission is expected to reduce this bound down to 2% [58J, 
if statistical anisotropy is not observed. Anisotropy in the non-Gaussianity can be dominant, 
which means that /nl may feature an angular modulation on the microwave sky. 

The above suggest that cosmological observations allow for detailed modelling and open a 
window to fundamental physics complementary to Earth based experiments such as the LHC. 
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